Astronomy & Astrophysics manuscript no. Fel44 


February 2, 2008 


(DOI: will be inserted by hand later) 





Detection of elements beyond the Ba-peak 
in VLT+UVES spectra of post-AGB stars* 

Maarten Reyniers and Hans Van Winckel** 



m 
O 
O 
(N 

< 



> 

00 

o 

oo 

o 
m 
o 



Oh' 
i 1 

O 1 



X 



Departement Natuurkunde en Sterrenkunde, K.U.Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium 
Received 14 May 2003 / Accepted 1 August 2003 

Abstract. In this letter, we report on our successful systematic search for lines of elements beyond the Ba-peak in spectra 
of s-process enriched post-AGB stars. Using newly released atomic data from both the VALD database and the D.R.EA.M. 
project, we could derive abundances for several elements heavier than europium for three objects, on the base of high quality 
VLT+UVES spectra. The abundances of these elements are of particular interest since they turn out to be powerful constraints 
for chemical evolutionary AGB models. Their high abundances indicate that, also in only moderately metal deficient AGB stars, 
production of lead is expected to be significant. 
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1. Introduction 

There is now general consensus that the neutron nucleosynthe- 
sis in AGB stars is triggered by prot on engulfment in the 12 C 
rich intershell (e.g. Bus so et alJl999t) . A significant 13 C-pocket 
is formed in a limited p/ 12 C re gime and the neutron irradiation 
occurs in radiative conditions JStraniero et alJl9 95) during the 
interpulse period, when this 13 C pocket is consumed by a cap- 
tures producing 16 and free neutrons. The dredge-ups asso- 
ciated with the thermal pulses during the AGB evolution will 
bring the enriched material to the stellar surface. There is, how- 
ever, no consensus on the mechanisms of both the dredge-up 
processes and the parti al mixing mechanism of protons in the 
intershell (e.g. Lug aro et all2003l and references therein). 

The characteristics of the s-process nucleosynthesis are ex- 
pected to be strongly dependent on the metalhcity since the 
12 C in the intershell is of primary origin, and thus largely in- 
dependent on the initial 12 C content. Moreover, since the p/ 12 C 
regi me for which the s-process is expected to be active, is lim- 
ited dGorielv & Mowlavill2000l) . the uncertainties on the neu- 
tron irradiation are mainly from the extend of the proton en- 
gulfed zone in the intershell, and much less on the ad-hoc as- 
sumptions of the proton profile. In reduced environments, there 
are more neutrons available per iron seed, pushing the s-process 
to heavier elements. The prediction that in metal deficient 
([Fe/H] < -1.0 ) stars, the productio n of Pb and B i is expected to 
be significant JGallino et alJll998t iGorielv & MowlavH EoOO) 
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was recently confirm ed by the detection of such Tead-stars' in 
extrinsic CH-stars bv lVan Eck et alJ fcOOll) . 

The simple picture, in which the metalhcity is the main 
parameter determining the s-process characteristics is, how- 
ever, much less clear after the finding of a number of objects 
of the same metalhcity, but with a large variety of s-process 
characteristics: s-process enriched objects of low metalhcity 
(-3 < [Fe/ H1<-1) were found without the exp ected Pb en- 
richment jAokiet alJl200lt I Van Eck etai1l2003l) : and in less 
deficient post-AGB objects, large overabundances were found 
in which again a large spread in s-process efficiency was ob- 
served at similar m etalhcity (I Van Winckel & Revniersl feoOO: 
iRevniers et all2003l) . 

In this paper we focus on a systematic search of the 
very heavy elements beyond the Ba-peak for the intrin- 
sic post-AGB stars. The reason for this search is twofold: 
an accurate distribution of s-process overabundances is a 
prime tool to constrain the s-process nucleosynthesis and, 
second, in two, somewhat hotter post-AGB stars, hafnium 
(Z = 72) lines were tentati vely detected, yielding a rela - 
tively large overabundance JVan Winckel & Revniersl 2000). 
Moreover, new atomic data for s-process elements beyond the 
Ba-peak b ecame recently ava ilable, either through the VALD 
database (Ku pkaetalJ Il999h . or via the new D.R.E.A.M. 
project (http://www.umh.ac.be/~astro/dream.shtml). We used 
our high-resolution, high signal-to-noise (S/N) VLT+UVES 
sp ectra for this purpose . These spectra are described in detail 
in Reynie rs et alJ l|2003i). Here, we only note that they cover a 
spectral interval from 3745 A to 10550 A; that their resolving 
power varies between 55,000 and 60,000; and that the S/N is 
lying between 100 and 200 (depending on wavelength). 
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2. Systematic search 

2.1. IRAS 06530-021 3 

We first concentrated on the spectrum of th e heavily en- 
riched IRAS 06530-0213 jRevniers et all2003l) . We extracted 
all atomic data for lines of elements between gadolinium 
(Gd, Z = 64) and platinum (Pt, Z = 78) from either the VALD 
database or the D.R.E.A.M. project. By assuming an ad-hoc 
high abundance, we determined a list of the strongest lines of 
all elements using a model atmosphere of IRAS 06530-0213. 
The eventual presence of those lines was checked on the spec- 
trum until the estimated strengths became too small to detect. 

After that, the detected lines were thoroughly checked in 
four different ways: (1) Lines with a suspected profile (asym- 
metry, continuum placement, cosmic hit, end of CCD) were im- 
mediately discarded; (2) Lines were compared to the spectrum 
of HR4656, a fast rotating (Vsinz*^ 210kms _1 ) Beta Cephei 
with spectral type B2IV, in order to check f or telluric blending . 
Also the solar spectrum identification bv iMoore et al.l £l966) 
was used for this purpose; (3) We also checked for possible 
diffuse interstellar bands (DIBs) in proximity of the line under 
study; (4) Blends are not always "visible" by a simple inspec- 
tion of the line profile. Therefore, a final check was performed 
on the lines by the use of the VALD database. We extracted all 
lines from VALD in the immediate vicinity (+1 A) of the line 
under study, and calculated the equivalent widths of all these 
extracted lines, using the abundances of IRAS 06530-0213. 
Such an exercise is highly necessary since it reveals sometimes 
strong blends (mostly from other s-process lines) which were 
not detected at first glance. 

The final selection for IRAS 06530-0213 contains three 
lines of singly ionised gadolinium (Gd n, Z = 64), two lines of 
singly ionised ytterbium (Ybn, Z = 70) and two surprisingly 
strong lines of singly ionised lutetium (Lun, Z = 71). We also 
found one line of tungsten (W, Z = 74), but the line is blended 
with a Sm n line. Nevertheless, we decided to derive a W abun- 
dance from this line by using spectrum synthesis (see Sect. [5). 
These lines are all shown in Fig.Q their atomic data and equiv- 
alent widths can be found in Tab.^ 

2.2. IRAS 05341 +0852 

All lines detected in IRAS 06530-0213 were obviously also 
detected in the VLT+UVES spectrum of IRAS 05341+0852, 
but, as can be seen on Fig.[0and Tab.[Q not all lines were used 
for an abundance determination since some are too blended 
with other s-process lines (e.g. the Gdn line at 5092.249 A) 
or some are even too strong (e.g. the blended Wn line at 
5 104.432 A). On the other hand, three additional Gdn lines 
and one additional Yb n line were found suitable for an abun- 
dance determination. Plots and atomic data can again be found 
in Fig.^and Tab. ^respectively. 

2.3. IRAS 081 43-4406 

The line strengths of the trans-Ba elements detected in 
IRAS 06530-0213 and IRAS 05341+0852 are much smaller 
for IRAS 8143-44 06, due to the less extreme enhancement of 
this object (Re vniers et al.l2003h . Some lines are even untrace- 



able (e.g. the Gdn line at 5469.713 A). As a result, only three 
lines were found suitable for an abundance determination. 

3. Abundances 

Tab. [2 also contains the abundances derived from the se- 
lected lines. These ab undances were calculated using Kurucz' 
model atmospheres kurucj 1 19931) and the latest version 
(April 2002) of Sneden's Stellar Line Analysis Program 
MOOG. Model parameters were taken from our earlier 
studies (IRAS 05341+0852: IVan Winckel & Revniers! l2000t 
IRAS 06530-0213 and IRAS 08143-4406: iRevniers et alJ 
l2003h . 

The W line at 5104.432 A does not fulfill the criteria men- 
tioned above since it is blended with a Sm n line at 5 104.479 A. 
In the case of IRAS 06530-0213, the total equivalent width 
of the blend is 73.6 mA, and the predicted contribution of the 
Sm ii line is 22 mA (assuming log e(Sm) = 2.09). Spectrum syn- 
thesi s with a m acroturbulent broadening £ m = 19.5 km s (see 
IRevniers et al.l2003l) yields a W abundance of log e(W) = 3.04. 
Since it concerns an isolated blend, we deduced the abundance 
also without profile match solely on the equivalent width ba- 
sis. This technique yields an abundance of log e(W) = 2.92. We 
prefer to use this lower abundance since it is independent of the 
adopted £ m . However, one has to treat this abundance with cau- 
tion since there might also be a blend present in the blue wing 
of the profile. We could not identify this blend, but it seems also 
present in the spectrum of the two other objects (see Fig. 
In the case of IRAS 05341+0852, we did not derive a W abun- 
dance from the blended line at 5 104.432 A, due to the relatively 
high contribution of the Sm n line in the blend (80.5 mA assum- 
ing log e(Sm) = 2.20). Also for IRAS 08143-4406, the relative 
contribution is considered to be too high for a reliable abun- 
dance determination. 

For Lu, hyperfine splitting (hfs) of the energy levels is 
expected. The Lu abundance was therefore deduced by spec- 
trum synthesis. For the line at 622 1.890 A, there are nine hfs 
components . Hfs A and B constants for the lower level were 
taken from iBrix & Kopfe rmannl lll952l) : for the upper level 
from lDen Hartog et alJ Jl998h . For the 6463. 107 A transition, 
there are only three components, and the splitting is entirely 
due to the lower level, since for the uppe r level 7 = 0. Hfs 
consta nts for this level were also taken from lDen Hartog et alJ 
(1998). Relati ve strengths were c alculated with the equations 
given in ICondon & Shortly! J 1 9 3 5b - The effect of hfs on the 
Lu abundance determination is considerable. In the case of 
IRAS 06530-0213, the difference in abundance between a hfs 
and a non-hfs synthesis (Ahfs) is -0.27 dex for the 6221.890 A 
line, and -0.14dex for the 6463. 107 A line. The effect is 
even larger for the cooler IRAS 05341+0852: Ahfs = -1.37 and 
-0.48 for the two lines respectively. For the hotter and less en- 
hanced IRAS 08143-4406, hfs only affects the line profile, but 
not the equivalent width. 

In Tab.|3]a summary of the final abundances is given. In or- 
der to compare these abundances with the solar ones, it should 
be checked that the solar abundances are based on the same 
sources of gf values as used in the present paper. This is, how- 
ever, for several reasons not always straightforward. The orig- 
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Table 1. Atomic data and abundances of individual lines of s-process elements with Z > 63. 
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IRAS 05341 


IRAS 06530 




IRAS 08143 


A 


X 


logs/ 


ref. 


W,i 


loge 




log 6 




loge 


(A) 


(eV) 






(mA) 




(mA) 




(mA) 




Gdn 


5092.249 


1.727 


-0.529 


(1) 






45.6 


2.72 


9.5 


2.04 


5108.903 


1.659 


-0.322 


(1) 


66.2 


2.26 


40.5 


2.39 


9.7 


1.78 


5357.794 


1.157 


-1.470 


(1) 


36.5 


2.59 










5469.713 


1.102 


-1.540 


(1) 


33.8 


2.56 










5733.852 


1.372 


-0.893 


(1) 


63.6 


2.52 


36.4 


2.64 






6080.641 


1.727 


-0.926 


(1) 


29.9 


2.42 










Ybn 


5352.954 


3.747 


-0.340 


(2) 






54.8 


2.59 






5651.988 


3.747 


-0.800 


(2) 


62.9 


2.88 


38.3 


2.84 






6489.298 


4.153 


-0.810 


(2) 


48.2 


3.07 










Lun 


6221.890 


1.542 


-0.760 


(3) 


245. 


1.58* 


187.4 


1.74* 


42.1 


0.97 


6463.107 


1.463 


-1.050 


(3) 


230. 


1.53* 


157.9 


1.78* 






Wn 


5104.432 


2.355 


-0.910 


(4) 






73.6 


2.92 1 " 







Refs. (1) Meggers et al. 1975 via Kurucz 1993 and VALD (2) Biemont et al. 1998 via D.R.E.A.M. 
(3) Ouinet et al. 1999 via D.R.E.A.M. (4) Clawson & Miller 1973 via Kurucz 1993 and VALD 
*abundance derived with spectrum synthesis (see text) 1 Sm blend taken into account (see text) 



inal re ference for the solar Gd abundance is Bergs trom et alJ 
dl988h and it is based on eight Gdn lines, without overlap 
with the lines in this paper. Ho wever, the log gf va lues of 
these eight lines are in our source (Meggers et al. 1975) some- 
what lower (on average 0.14dex) than the Bergstrom et alJ val- 
ues. This might be an indication that me lMeggerTetalT val- 
ues are systematically too low and that the Gd abundances 
should be lowered by ~0.14dex to be consistent with the so- 
lar value. For Yb, the details of the solar abundance calcu- 
lation are unfortunately lost (Grevesse, private communica- 
tion). ^_TJie_s^>hu~ Luabjmda^^ the 6221 A 
line dBord et alJll998t ben Hartog et al.lll998l) . with the same 
log gf and hfs parameters as in this paper. Finally, the situation 
is less clear for W, since the solar value is derived from neu - 
tral W lines dHolweeer & Wernerll982llDuauette et al.ll98ll) . 
Moreover, there is a large discrepancy (0.42 dex) between the 
photospheric and the meteoritic value fo r this element, the 
photo spheric being the most uncertain one JGrevesse & Sauval 



4. Discussion 

To our knowledge, this is the first systematic search for trans- 
Ba elements in intrinsically enriched objects. Especially the 
detection of two strong Lu lines is astonishing, but an alter- 
native identification was not found. Tomkin & Lambert ( 1983) 
reported the Lun 622 1.890 A line to be blended in the clas- 
sical Ba-star HR774, but a blending species was, however, not 
given. We think this detection is real, and that the non-detection 
by previous authors is in fact a consequence of the lack of suit- 
able line lists in the past, the time consuming procedure of a 
systematic approach in this matter and the lack of high S/N 



Table 2. Final abundances of elements beyond the Ba- 
peak. For each ion, we list the number of lines N, 
the mean equivalent width Wx, the absolute abundance 
loge (i.e. relative to H: log e = logX/H + 12), the line-to- 
line scatter cr and the abundance relative to iron [el/Fe]. 
ss stands for spectrum synthesis, bl means corrected for 
blend. (Photosph eric) solar abundance s are taken from 
the review bv iGrevesse & Sauvall 415)98): log e(Gd) = 1.12, 
log e(Yb) = 1 .08, log e(Lu) = 0.06, log e(W) = 1.11. 



ion 


N 


W A log e cr 




[el/Fe] 






IRAS 05341+0852 






[Fe/H] 




-0.7, T cff = 6500K, lo£ 


\g 


= 1.0 


Gdn 


5 


46 2.47 0.13 




2.07 


Ybn 


2 


56 2.98 0.13 




2.62 


Lun 


2 


ss 1.56 0.04 




2.22 






IRAS 06530-0213 






[Fe/H] 




-0.5, T cff = 7250 K, logg 


= 1.0 


Gdn 


3 


41 2.58 0.17 




1.92 


Ybn 


2 


47 2.71 0.18 




2.09 


Lun 


2 


ss 1.76 0.03 




2.16 


Wn 


1 


bl 2.92 




2.27 






IRAS 08143-4406 






[Fe/H] 




-0.4, T cff = 7250K, lo£ 


\g 


= 1.5 


Gdn 


2 


10 1.91 0.18 




1.18 


Lun 


1 


42 0.97 




1.30 



spectra obtained with efficient spectrographs on large aperture 
telescopes. 

As a final check for the identifications presented above, 
we compare the derived abundances with the predictions of 
the chemical evolutionary AGB models of the Torino group 
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Fig. 1. Detection of lines from elements beyond the Ba-peak in the spectra of IRAS 05341+0852, IRAS 06530-0213 and 
IRAS 08143-4406. All these lines are intensively checked for blends, telluric pollution and DIB contamination (see text for de- 
tails). From top to bottom are plotted: IRAS 05341+0852, IRAS 06530-0213 and IRAS 08143-4406. The spectral lines marked 
with an arrow are used for abundance determination. 



(e.g. lGallino et alll998llBusso et all 199911200 lb . These mod- 
els are frequently used to fit observed abundances by adjusting 
the 13 C pocket till the predicted and observed s-process dis- 
tributions match. The fitting procedure, together with a more 
ge neral discussion of th e comparison, is described in detail 
in R evniers et alJ d2003h . Here, we focus on the elements be- 
yond the Ba-peak. In Fig. |2]we graphically compare the de- 
rived abundances with the best-fit model. The parameters of 
each model are given in the upper left corner; the initial mass 
is 1 .5 M for each model. 

The predictions for the elements beyond the Ba-peak are 
consistent with the observed abundances. This is an additional 
argument strengthening the line identification for these ele- 
ments. Only the Yb abundance of IRAS 05341+0852 is quite 
off the model curve (by +0.3 dex). More specifically, there 
is a large difference between the Yb abundance and the Lu 
abundance for IRAS 05341+0852, which is not predicted by 
the models. This difference is not seen in the abundances of 
IRAS 06530-0213, although three lines out of five are in com- 
mon for these elements (see Tab.^. 

As can be seen from the comparison between two models 
with a different I3 C pocket for IRAS 06530-0213, the abun- 



dances of the very heavy s-process elements are most sensitive 
to the adopted I3 C pocket, Hence, these elements (especially 
Hf and W) are ideally suited to discriminate between possible 
13 C pockets. Unfortunately, a Hf abundance is very difficult to 
obtain for these objects since it has only suitable lines in the 
blue, and the W abundance derived from the 5 104.432 A line 
depends on the adopted Sm abundance (see Sect. [3}. 

We conclude that we found trans-Ba abundances in strongly 
enriched post-AGB stars. Despite their difficult abundance de- 
termination, they are sensitive indicators for the strength of the 
neutron irradiation during the AGB interpulse phase. Note that 
also in these post-AGB stars with large overabundances of s- 
process elements beyond the barium stability peak, a strong 
overabundance of lead (Pb, Z = 82) is expected. Unfortunately 
Pb has no suitable non-blended lines in the optical spectrum 
for these stars. Our detection of very heavy s-process elements 
indicate that also for moderate metallicities, Pb production is 
expected to be significant. 
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The anonymous referee is thanked for the suggestion of the impor- 
tance of the Lu hyperfine spitting; Andrew McWilliam for his help 
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Fig. 2. Comparison of the observed abundances with the best 
fit Torino model predictions. The parameters of the adopted 
model are given in the upper left corner. Note that the spe- 
cific model is obtained by fitting the Zr and La abundances 
(see lRevniers et al.l2003l for details), so without special atten- 
tion to the very heavy s-process elements. The errorbars plot- 
ted on the figure are the line-to-line scatters that are given in 
Tab. |3 The actual error on an abundance is, however, in most 
cases considerably larger than this value, due to uncertainties in 
the atmospheric parameters, undetected blends, uncertain solar 
abundances, etc. . . One example is the solar abundance of W, 
for which there is a difference o f 0.42 dex between t he photo- 
spheric and the meteoritic value JGrevesse & Sauvall9 98). 
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